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instant or when it does not. Particular attention is paid to the corrections ensuing from the non- 
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classical Erlang distribution is considered in detail. 
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I. INTRODUCTION 

In this article we study exit times of continuous-time 
random walks (CTRWs) with drift. By this we under- 
stand a random process X t whose evolution in time can 
be thought of as the result of the combined effect of a con- 
stant drift and the occurrence of random jumps. Thus, 
for t 5s to wc define 



X t = X tn + v{t - t Q ) + S t , S t 



Tl = l 



Jn 6{t-t n ), (1) 



where 9(u) = 1 for u ^ 0, t n = t\ + • • • + T n are the jump 
times, J n the jump magnitude and v > by convention. 
Note that as the notation suggests, to coincides with a 
jump. The process X t may represent physically the posi- 
tion of a diffusing particle, for instance. All through this 
paper we will assume that (i) the sojourn times r„ > 
are independent and identically distributed (i.i.d.) ran- 
dom variables with probability density and cumulative 
distribution function (PDF and, respectively, CDF) ip(t) 
and $(t) = J Q ip(t')dt'; (ii) {J„} is a sequence of i.i.d. 
random variables with common PDF /i(-); and (hi), that 
J rn is independent of r„ for any n, m. 

In absence of drift, and when the holding times r n are 
exponentially distributed, t„ ~ £(A) for some A > 0, the 
jump process St is a classical compound Poisson process 
(CPP) characterized by having independent increments 
St+h — St; in addition the associated "arrival-process" 
N t = Y,n=i °( t ~ *«) is Poisson distributed: N t ~ V(Xt). 
Therefore, drift-less CTRWs can be considered as a gen- 
eralization of CPPs where the holding times of the pro- 
cesses St are arbitrarily distributed. In statistical physics 
such drift-less CTRWs have been widely used after the 



work of Montroll & Weiss [l|,[2| and their relevant statisti- 
cal properties, like correlation functions and the behavior 
in the continuum limit, 1 a subject of intense interest 0- 
@- Applications of CTRWs can be found in the study 
of transport in disordered med ia ( e.g., anomalous 
relaxation in polymer chains [10(, sandpile and earth- 
quake modeling (e.g., [U fl2|). random networks [Hj], 
self-organized criticality in gran ular systems fl4j . scal- 
ing properties of Levy wal ks Il5ll. electron tunneling [l6| , 
transmission to mog raphy [171 Il8j , distribution of matter 
in the universe [19( and changes in stock markets due to 
unexpected catastrophes [2(| . More recently, the use of 
CTRWs has been advocated to give a microscopic, tick- 
by-tick, description of financial markets: see [2ll - l25l ]. A 
comprehensive review of CTRW applications in finance 
and economics is given in (26j . 

Physically, the introduction of these general CTRWs 
stems from the fact that in many settings the exponential 
holding-time assumption may be inadequate to describe 
the physical situation — see [HI, [2l[ . Additional motiva- 
tion arises from the observation that there is an extensive 
number of different physical systems that show some sort 
of anomalous diffusion — a subject of great interest from 
the viewpoint of statistical physics in the last years — , 
and that this anomalous behavior can well described by 
fractional Fokkcr-Planck equations obtained after impos- 
ing the continuum limit on CTRW models (27l - l29| . 



The further addition of the drift term to a CTRW, as 
we do in Eq. ([1]), is a natural and significant incorpora- 
tion. The resulting process — which can be viewed as the 
discrete analogue of a (fractional) diffusion with drift — 
is known for playing a fundamental role in the modeling 
of the cash flow at an insurance company [30l l3l| and, 
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1 This limit corresponds to the assumption that the sojourn time 
between consecutive jumps and their characteristic size go to zero 
in an appropriate and coordinated way. We will show an explicit 
example of this limit in the Appendix. 
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more recently, it has been shown that it also rules the 
rate of energy dissipation in nonlinear optical fibers [32| . 
Indeed, present and forthcoming results can be also of 
interest in transport in amorphous media l33|, models of 
decision and response time in psychology |34| and neuron 
dynamics [35| . 

In all these scenarios one is faced with the basic prob- 
lem of determining the first-passage time for a CTRW 
with drift, a question that has been theoretically consid- 
ered in the past [36U431 ] . The usual approach taken there 
entails the computation of transition and first-passage 
time distributions of the process, e.g. [13, EH, an d the 
results are typically obtained under the assumption of 
the continuum limit [40l - [43| . This procedure has incon- 
testable advantages for obtaining the leading-order be- 
havior, but it is not adequate for a detailed analysis of 
the statistical properties of the process at the inter-jump 
timescale, as in [|1 HI . 

Motivated by the above, here we pose the problem of 
evaluating the mean exit time from the interval (0, b) of 
a drifting CTRW X t when the only available informa- 
tion is the present state X r = x, where r is the 
present chronological time — note that, by adjusting the 
time clock and spatial scale, the results carry over to any 
interval (a, b) and initial time to- Here X t is given by 
Eq. ^ where St is the jump part and N t the associated 
counting process (a renewal process). In the drift-less 
case previous work in this regard includes that of [22], HH 
where a linear integral equation for the mean escape time 
after a jump off a given interval is derived. We note 
however that these results do not cover a generic situ- 
ation. Indeed, while for CPPs (as for the general Levy 
processes) Markov property implies that results derived 
starting at a jump time carry over to arbitrary present, 
no such inference is possible for a generic CTRW due to 
its non-Markovian nature. Thus, escape times depend 
on the actual state and available information and hence 
the question as to how to generalize the former results 
to general present time r appears naturally. We remark 
that implications ensuing from the lack of Markovianess 
have been ignored at large in the literature, a gap that 
we have intended to fill in — see [44j|. In particular, it 
remains an open question to what extent dropping the 
assumption that "the present is a jump time" affects the 
relevant probabilities. Here we address these issues and 
generalize the results of [H, [23| in a twofold way by as- 
suming that (i) a drift vt > operates on the system 
and (ii) the present r is an arbitrary time, not necessar- 
ily a jump instant, and the observer has knowledge of the 
present, but not of the history, of the system. 

The interest of this problem goes far beyond the purely 
academic since such a situation may appear in several 
different physical contexts. For example, one might be 
interested in predicting the mean time for an insur- 
ance/financial company to go bankrupt from the knowl- 
edge of just the actual company budget, i.e., when the 
information regarding the company's past performance 
has not been disclosed. A second example is provided 



by the study of the distribution of inhomogeneitics in an 
optical fiber; it was found [32[ that the signal's energy 
amplitude at a point t (t is the spatial variable in this 
setting) involves a CTRW with drift X t ; in this context 2 
one typically knows only the value of the energy at the 
observation point r, not on the whole fiber. More gener- 
ally the approach will be relevant in situations where ei- 
ther the elapsed time between events is "large" (it might 
be as large as years, in a context of catastrophes obser- 
vation) or when the event's times t n are not physically 
measurable observables and only mean escape times are. 
(Note that both the mean exit time and initial time are 
typically macroscopic magnitudes.) 

The article is structured as follows. In Section |TT] we 
show how the solution to these problems involves ideas 
drawn from renewal theory and solve the simpler case 
when v = 0. The case v > is considered in the next sec- 
tions where it is found that key properties of the obtained 
equation depend on the sign of the jumps. In Section HTT1 
it is shown that if this sign is positive the solution can 
be given in closed form by Laplace transformation — cf. 
Eqs. (fT2|) and (fT3|). For the case when drift and jumps 
have opposite signs we find integral equations that the 
relevant objects satisfy, see Section ITVl but no closed so- 
lution can be given in a general situation. Section [V] 
addresses the most general scenario in which jumps J n 
can take both signs. Solvable cases are discussed there. 

In all cases we exemplify our results by considering 
the particular instance when sojourn times have Erlang 
distribution, £r{\, 2). £r{\, n) corresponds to having a 
sum of n independent exponential variables and hence 
generalizes the exponential density in a natural way, 

\nxn— 1 

(n- 1)! 

while it maintains an adequate capability to fit mea- 
sured data. From a physical perspective these facts 
make this density a natural candidate to describe multi- 
component systems which operate only when several 
independent, exponentially distributed operations have 
been completed or whenever there is a hidden Poissonian 
flux of information and jumps only appear as the out- 
come of two or more consecutive arrivals. This explains 
the interest that it has drawn in the field of information 
traffic [45], [4(| . Similarly the appearance of this distri- 
bution to model transaction orders in financial markets 
can also be expected since it takes, at least, two arrivals 
(buy and sell orders) for a transaction to be completed. 
For further applications to ruin problems and insurance 
see [13, S3. 



2 This situation where t represents the the space variable may 
occur in different physical systems, with X standing for some 
physical observable of interest, like the energy. In this connec- 
tion, r might represent the location of a detector, or a sensitive 
part of the appliance at which measures are taken. 
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II. THE PROBLEM 

Recall that we aim to study exit times of a drifting 
CTRW X t given the present state X r = x. To this end 
let r+tl' r be the first time past r at which X t exits (0, 6); 
T b (x,r) be its expected value: T b (x,r) = E[t^' r ]; finally 
let Tb{x) denote the mean exit time off (0, b) after a jump 
t n (loosely one has Tb(x) = T b {x, t n ), n = 0, . . . , oo; how- 
ever the relation between both quantities is not trivial, 
as we see below in Section [TO]) . Note also that here and 
elsewhere we use E[-] to denote expectation. 

In the exponential Markov case t^' r is independent 
of r, T b (x,r) = T b {x) and it only remains to formulate 
(and solve) the equation that this object satisfies. How- 
ever this situation no longer holds in the generic, non- 
Markovian case where t£ ' r does depend on r. We find (see 
below and Section HTl]) that the relation between T b {x,r) 
and Tf,(x) involves the "excess life" E r = ij\r r +i — r, or 
time elapsed until the next arrival occurs. We now sketch 
classical renewal theory (see [4!|[5(j and [j|) that shows 
how to construct the CDF $(t|r) = ¥{E r ^ t} of E r , 

Let m(t) = E[7V t ] be the mean number of jumps up to 
t: the renewal function. It satisfies the integral renewal 
equation 



m(t) = #(t) + / m(t - t')ip(t')dt' . 



(3) 



Then, by using the total probability theorem it can be 
proved that 



$(i|r) 



r+t 



[l-V(r + t-t')}dm(t'). (4) 



Upon solution of the above integral equations we obtain 
<£(i|r). Actually, they can be solved with all generality 
by recourse to Laplace transformation. Let g(s) be the 
Laplace transform of a function g(t) so that 



1 

3 2ni 7 C _ ?; 



c+ioo 



e st g(s)ds, c > 0. 



(5) 



Then, Eqs. © and (TJ| allow to recover the distribution 
of E r in closed form via 



1 V(s) 



l-V(s)' 



(6a) 



l-tp(s) / e- sl m{l)ldl, (6b) 



rh(s) 



</>(s\r) 



where rh(s) and </>(s|r) are the Laplace transforms of m(t) 
and <f>(t\r) = d t $(t\r). 

If v = these expressions can be used to relate Tb(x, r) 
and Tb{x). Indeed let t^ r = t n , say, be the "last" jump 
time and E~ = r — t n the time elapsed from t n to the 
present. Then, with £"+ = E r one obviously has that 
E~ + £"+ = t n+ i — t n = and that the exit time 

right after t n is that after r, t^' r , plus E~ , and hence 



T b (x) = e[e; 



Eft; 



= E[T n+ i]-E[ J E7+]+T 6 (a;,r) 



Thus, Tb(x,r) follows adding a correction term to Tf,(x) 
which depends only on r. Finally Tb(x) is obtained by 
solving a linear integral equation — see [H, [23|. Un- 
fortunately when v ^ this simple argument fails as 
then knowledge of the present position does not entail its 
knowledge at t„ . In the next sections we derive the rel- 
evant correction to the mean exit time. This correction 
depends now in all parameters r, x and b — see Eqs. © 
and (|26p below. 



III. JUMP PROCESS WITH FAVORABLE 
DRIFT 

In this section we consider the case when both drift and 
jumps have a positive sign, i.e. when X t is increasing. As 
a result, the process can only leave the interval through 
the upper boundary b. Let us assume that at time t = r 
the system is in X r = x € (0,6), and that the excess 
life E r is known in advance, E r = I. If the excess life is 
longer than g = < I, the drift will drive the process 
out of the region at time r + g, before the next jump 
takes place at tjy r +i. Conversely, if I ^ g at least a 
jump of size Jn t +i = u will occur prior to exiting the 
interval. Note that just before tN r +i the process X t is 
no longer at x, but at x + vl. Therefore, two possible 
scenarios appear: either the jump size is larger than the 
remaining distance up to the upper boundary, u b — 
x — vl, and the process leaves the interval at r + E r , or 
it does not. In the latter case the problem renews from 
t = r + I, X r +i = x + vl + u, so the mean escape time 
will be increased by an amount T& (x + vl + u). It can 
be proven that these considerations imply that Tb(x,r) 
must satisfy — recall that tp(-), h(-) are the waiting-time 
and, respectively, jump PDFs — the following equation 

re roo 

T b (x, r) = [1 - $ (g\r)} g+ lcf)(l\r) / h(u)dudl 

JO Jb-x-vl 

r-h — X — vl 

4>{l\r) / h(u)[l + T b (x + vl + u)]dudl,(8) 

'o Jo 

an expression that relates Tb(x,r) and Tb(x). Note that 
since $(Z|r) depends on r so it does the mean time 
Tb(x,r). Finally, after some rearrangement, Eq. © can 
be conveniently written as 



Tb{x,r) 



[l-$(l\r)]dl 



+ - / (f>\Q--r) h(u)T b (b- z + u)dudz. (9) 

V In V V J ,/n 



Similarly, by letting r — > 0, we find that T b (x) must 
satisfy 



T & (x) 



[1 - *(Z)] dl 



(7) 



i r e f z\ r z 

+ - il>\Q--) h(u)T b {b- z + u) 
v Jo v vJ J 



dudz, (10) 
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X G (0,6). Equation © along with (|10l) allow us to 
solve the posed problem. The second of these defines an 
integral equation for T b {x) which, if v = 0, reduces upon 
appropriate change to that of [U [23| . It is remarkable 
that Eq. (|10p can 6e solved in a fully explicit way. To 
this end we define the allied object F(y), as the solution 
of the following integral equation 



F(y) 



y/v 



y - z 



[1 - *(J)]tB 



/i(u)F (z - u) dudz, (11) 



for y G R+. Then it follows that T b (x) = F(6 - x), 
for x G (0,6). We note further that taking a Laplace 
transformation in Eq. pip we find that 



F(s) 



1 



1 — ip(sv) 



vs 2 l _ ip(sv)h(s) 



(12) 



where F(s) = J °° F(y)e-'*dy, = J °° /i(u) e -™du, 
and ?/>(s) = J °° tp(s)e~ st dt. Here s = sji + isi is complex 
and Sfl > 0. Further, we also have 



J(s\r) =F(s) 



1 - ft(s) 



1 — ip(sv)h(s] 



4>{sv\r) — ip(sv) 



(13) 

where the correction to the mean time after a jump is 
clearly displayed and again, for convenience, we defined 
J(y\r) = T b (b - y\r), and extended © to y ^ 0. Thus 
F(y) and J(y|r) can be recovered by Laplace inversion, 
cf. Eq. ([5]), and Tf,(x) and T b (x\r) will eventually follow. 

Equation (|13p has several mathematical limits of ap- 
plied interest. We hrst consider the situation when v is 
small. Let fi = E[r„ + i], /i r = E[_E. r ]; then, using that 
ijj(sv) = 1 — fisv + 0(v 2 ) and so forth we see that for 
small v, J(s\r) has an expansion in powers of v as 



J(s\r) = F(s) - 



Mf) 

1 - h(s 



■(fj, - fi r )fJ,V + 0(v 2 ), 

(14) 

which implies, in particular, Eq. ([?])■ Another interesting 
case is obtained letting r — > oo: the steady-state solu- 
tion. This limit is relevant since it can be associated 
to a situation in which the only information available to 
the observer is the present value of the stochastic process, 
not even the starting point. Recalling that by the renewal 
theorem lim m(t)/t = l//i, Eqs. (jB]) and (fT3"|) yield 

t— ¥oo 



</>(s|oo) 



i-#0 



SjJL 



(15) 



Moo) = F ( ,s)-i#i-ii±^^.(16) 

We illustrate the ideas above by detailing the case 
when jump magnitudes have an exponential distribution 



h(u) = je 7 ", where 7 > is a real parameter, and so- 
journ times an Erlang distribution r„ ~ £r(A, 2). Hence 



h(s) 



7 



7 + s 



X 2 



(17) 



In this case Eq. (fTTj) yields that F(s) is the following 
rational function 



F(s) 



1 



2A7W + v(2X + -yv)s + v 2 s 2 



vs 2 A(A + 2jv) + v(2X + jv)s + v 2 s 2 
Hence, the mean escape time of (0, 6) after a jump is 

27V 



(18) 



A 2 



A + 2 7 , 
l_ e -*-e i_ e -^+e 



Q 



(19) 



where we recall that vg = b — x is the initial distance to 
the boundary 6 and 



(20) 



with Re[z±] > 0. 

The evaluation of the mean exit time starting at r, 
T b (x, r), involves the renewal function and excess-life dis- 
tribution. We hrst obtain from Eq. ([6]) 




m(s) = 



A 2 



s 2 (2A + ,s) : 



(s\r) =tp(s) + 



(1 - e- 2Xr )s\ 
2{\ + s) 



By inversion we get m(t) = (2\t + e 2At — l) /4 and 



(21) 



$(i|r) = 1 - e 



-A/ 



1 + e 



-2Ar 



Af 



(22) 



Then by using Eq. (fl~3f we find at last that 



A 1 - e 



-2Ar 



2 z + — z_ 



T b (x,r) = T 6 (x) 
l_ e -*-e i_ e -*+e" 



(23) 



Plots of this function in terms of x are given in Fig. [T] 
for several values of r and a certain choice of the rest of 
parameters. 

Finally note how, in particular, if v = then 

T 6 (x) = I [1 + 7(6 - x)] = E[r„] + ,(24) 

1 - e~ 2Ar 
T 6 (x,r) = T 6 (x) . (25) 



IV. THE CASE OF OPPOSITE DRIFT AND 
JUMPS 

We now consider the case when the sign of the drift 
vt > is opposite to that of the jumps. In this case the 
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Figure 1: (Color online) Sample representations of T b (x,r) 
for different values of r: r = 0, i.e. Tb(x), solid (black) line, 
r = 0.4, dashed (red) line and r — 10.0, dot-dashed (blue) 
line. The rest of the parameters were chosen as follows: b — 
1.0, v = 0.1, A = 1.0 and 7 = 0.1. 



process can leave the interval through both the upper 
and the lower boundaries: the drift pushes steadily the 
process up, whereas the jumps threaten the system with 
a downside exit. The resulting process is a prototype 
model in risk management to describe the dynamics of 
the cahsflow X t at an insurance company under the as- 
sumption that premiums are received at a constant rate 
v > and that the company incurs in losses J n < from 
claims reported at times t n ,n = l,...oo (the Cramer- 
Lundberg model). 

As before, we analyze the evolution starting at t = r 
with X r = x. Then if E r = I, I > g, the drift will drive 
the process out of the region through the upper bound- 
ary at time r + g. Otherwise (I g) at least a jump 
Jn t +i = —% say, occurs at time ij\r r +i before escape, and 
two possible scenarios appear depending on the relative 
magnitudes of the jump size and the location of the pro- 
cess right before the jump, X t = x + vl: for u ^ x + vl the 
process will leave the interval through the lower boundary 
at r + E r , but when u < x + vl the process after the jump 
will remain inside the interval, X r+ i = x + vl — u, the 
exit problem will start afresh, and the mean escape time 
will be increased by T b (x + vl — u). Again these consid- 
erations imply that T b (x,r) and Tf,(x) must satisfy for 
< x < b 



T b (x,r) = [l-$(g\r)]g + l(f>(l\r) / h{u)dudl 

JO J x+vl 

PQ rX + vl 

+ (f>(l\r) h(u)[l + T b (x + vl-u)}dudl 

Jo Jo 

re 

[l-$(l\r)]dl 

v .L \ v 



r ) j h(u)Tb (z — u) dudz, 







T 6 (a;)= / [1 -*(*)]<« 
Jo 



h(u)T b (z — u) dudz. 



(27) 

Hence T b (x,r) follows in terms of quadratures also 
in this case, given Tj,(x). Unfortunately, unlike what 
happens for the case considered in the previous Section, 
Eq. (|27[) can not be solved in closed form for arbitrary 
PDFs tp(-) and h(-). Further progress can be made for 
Erlang times, £r(X, 2). Indeed, in this case (|27p reads 

2 - (2 + Xg) e~ Ae 



T b (x) 



A 



x ^ e -M,-x)/v / h{u)T b (z-u) dudz, (28) 



A 2 f\ 
"7a / ( z 



and hence, by differentiation we find that T b (x), for 
< x < b, also satisfies the following integral-differential 
equation: 

2A 



0\ \2 

T b \x) - -T' b (x) + -T b (x) 



A 2 

H — 2 / h(u)T b (x — u) du, 



(29) 



(30) 



subject to the following boundary conditions: 

lim Tft(x) = 0, and lim T b (x) = —1/v. 

x^b x^b 

We first consider a general solution to this equation 
extending it to the full real axis, so we will drop the 
subscript in T b (x). We find a solution by Laplace trans- 
formation as 

f 2A/, + (B,-2AA) P + **, 

A 2 [l -h(s)} -2\vs + v 2 s 2 

where A and B are T(0) and T'(0) respectively. By in- 
version, cf. Eq. ((5]), T(x) follows in terms of A and B. By 
requiring (|30[) we obtain a linear algebraic system for A 
and _B, which upon solution yields T b (x) in closed form. 

To be specific we consider the case when jumps are 
also exponentially distributed: h(x) = r ye~ 1 ' x . Then we 
have (|17p and T(s) is the rational function 



T(s) = 



(7 + s)[2A + (Bv - 2XA)vs + Av 2 
s 2 [A(A - 2'fv) + (<yv - 2X)vs + v 2 



(32) 



Upon re-scale of constants the inverse Laplace transform 
of (El) reads 



A 



A + 7X 



A 



B 



^(X-XA + vB)--^ 



(X-XA + vB) - f £ 



— 7«J T(x) 

e t+x _ I 

e+(e+-e-) 

et~ x - 1 



(26) 



(33) 
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with 



s v 2 2 V jv 



(34) 



Unfortunately the final expressions for A and B af- 
ter imposing (|30[) are not very illuminating so we do not 
transcribe them here. Sampling values for different pa- 
rameter specifications can be found in Fig. [21 

The limit b — > oo is interesting as T oc (x) gives the 
probability that X t ever hits 0. This corresponds to the 
classical ruin probability in an insurance context. It turns 
out that Too (a;) can be determined in a direct way that 
avoids solving the aforementioned linear system. With- 
out proof 3 we note that if A > 2jv then 



Too (^) 



2(1 + 

A - 2jv 



(1 + 72:) 



1 



E[ Tn ] E[J„ 

(35) 

while Too (a:) = oo otherwise. Once Tf,(x) is known, 
Tb{x,t) follows again by integration — see Fig. [2j 

V. TWO-SIDED JUMP PROCESS WITH DRIFT 

We finally consider the general case corresponding to 
a jump process where J n can take both signs and hence 
can exit (0, b) through either of the boundaries. The 
relevant analysis is similar to that of the last section if 
one incorporates the possibility of an upper exit due to a 
jump: If E r = I > g the drift drives directly the process 
through the upper boundary. If I ^ g at least a jump 
occurs before the exit of the process. When the size u of 
this jump is positive and u b — x — vl the process leaves 
the interval at time r + E r through the upper boundary, 
when u ^ — (x + vl) < the exit takes place through 
the lower boundary, otherwise the process remains inside 
the interval (0, b) and the problem restarts. Skipping 
minor details we obtain that T&(x) and T b (x,r) satisfy 
respectively 



T b (x,r) = [ B [l-<f>(l\r)}dl 
Jo 



b-z 



h(u)T b (z + u) dudz, (36) 



and 



T b {x) = j e [l-^(V)]dl 
Jo 



-M— 



b-z 



h(u)T b (z + u) dudz. (37) 



In a general situation, the latter integral equation is not 
solvable in closed form. To gain some insight we use the 



(a) 




3 We elaborate on a similar expression in the next section. 



(b) 



Figure 2: (Color online) Sample representations of Tb(x,r) 
for different values of r: r = 0, i.e. Tt(x), solid (black) 
line, r — 0.4, dashed (red) line and r = 10.0, dot-dashed 
(blue) line. The rest of the parameters were chosen as fol- 
lows, b = 1.0, v = 0.1, A = 1.0 and (a) 7 = 0.1; (b) 7 = 4.0. 
We see how, when drift and jumps have a different sign, an 
interesting behavior may appear in Tb(x,r). In the upper 
panel we observe that this function is no longer decreasing 
with x and has a maximum in the interior of the interval, cf. 
Fig. [T] In the lower panel we observe how a cross-over phe- 
nomenon may eventually take place reflecting the fact that 
Tb(x) need not be greater than Tb(x,r): this behavior can 
be traced to the fact that as r increases the jump probability 
increases, which results in a reduction of the escape time if 
the process is in the vicinity of the lower boundary but in an 
increment when close to the upper boundary. The maxima 
position moves toward greater values of x for larger choices of 
r, in both cases. 



decomposition h(x) = qh + (x) + ph-(x) where p is the 
probability that a given jump be negative, q = 1 — p and 
h±(-) are the jump PDF in the positive/negative regions, 
i.e. h+(x) = h(x)6(x) I q, say. 

Note that if h-(x) = 0,Vx ^ 0, we recover the case 
considered in Section Mil solvable via Laplace transfor- 
mation. It turns out that we can still construct an ana- 
lytic closed solution in the more general case when h- (x) 
vanishes only on (—6, 0) — but not on (—00, —b]. In such 
a case a negative jump will drive the process out of the 
interval (0, b) through the lower boundary. Thus p is re- 
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lated to the ruin risk in a economic scenario or to the 
breakdown probability in a physical system. The equa- 
tion for Tfc(x) reads in this case 



q - / i>( — 

v J- V v 



b- 



(x) = r [i - * (/)] di 

Jo 

h+(u)Tb (z + u) dudz, (38) 



while a similar expression, with <fr(-\r) in place of 
holds for Tb(x, r). Note that these equations are indepen- 
dent of the form of h-(-) and apart from the factor q in 
front of the second term they resemble Eqs. ([5]) and (ITU1) ; 
it follows that we can resort to the same technique used 
in Section irm We consider 



F(y) 



y/v 



v Jo V v 



[i - m{i)]di 



h(u)W (z - u) dudz, (39) 



for y G M. + which is again solvable by means of a Laplace 
transform; then T&(x) follows from Tb(x) = F(£> — x) for 
x G (0,6). With the previous selection for h(-) and ?/>(•), 
Eq. (|17[) . and in terms of z = sw we find 



F(s) 



v 2A7W + (2A + 7^)2 + z 2 

z p\ 2 -fv + A(A + 27«)z + (2A + -fv)z 2 + z c 



.(40) 



We first consider the case when A = jv. Under this as- 
sumption F(s) has poles at z — Zj = A (^q 1 / a e 2vij / 3 - l), 
j G {1,2,3}, and z = 0. Inverse Laplace transformation 
yields the mean exit time as 



1 -)- q-l/3 e -2nij/3 
( 

3zj 



(41) 



Returning to the general case we see that the inver- 
sion of the Laplace transform involves solving a cubic 
equation, and though explicit formulas are available the 
resulting expression is awkward. Still, the large b limit 
can be discerned with all generality. To this end we note 
that by appealing to Hurwitz's stability criteria it can 
be proven that all three roots Zj, j G {1,2,3}, of the 
denominator in expression f|40[) — apart from z = — 
have negative real parts. Hence, evaluating the inversion 
integral by residues we find 



T 6 (x) 



2 

PA 



(42) 



where Cj are certain constants. Thus, letting b — > 00 we 

see that T 6 (x) ^ ^ = pjfe^ ■ 

The evaluation of the correction to Tf,(x,r) does not 
present particular difficulties. We leave it as an exercise 
to the interested reader. 



VI. CONCLUSIONS 

We have analyzed the mean exit time for a general 
CTRW with drift. If the present coincides with a jump 
time we find that it satisfies a certain integral equa- 
tion whose solvability is analyzed. We consider next the 
generic case when the present is an arbitrary instant and 
the history of the system is not available to the observer 
and only the present state is. It turns out that the cor- 
responding escape time can be obtained by incorporat- 
ing an appropriate correction, which can be described in 
terms of the "excess life", a familiar object in renewal 
theory. We find that when the drift and jump compo- 
nents have the same sign the equations that these objects 
satisfy can be solved in closed form via Laplace transfor- 
mation, irrespective of the distribution; otherwise, one 
must restrict to particular choices of the sojourn-time 
distribution. The case corresponding to the classical Er- 
lang distribution is analyzed in detail. The more gen- 
eral case when jumps take both positive and negative 
signs is also considered and solved under certain severe 
conditions. We plan to generalize these ideas to a more 
general class of waiting-time distributions and pinpoint 
conditions that guarantee the reducibility of the original 
formulation to simpler differential equations. 

The relevance of these results from a physical perspec- 
tive is discussed in several connections of interest includ- 
ing possible applications to risk, finance and distribution 
of energy in optical systems, which will be the matter of 
future publications. We also point out the relevance of 
the approach whenever the time between events is "large" 
or when the arrival times are not physical observables. 
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Appendix A: Some remarks on the continuum limit 
and its relationship with fractional diffusions 

In this appendix we sketch how the approach followed 
in this article, based in the use of renewal theory, com- 
pares with the most traditional one which relies on the 
previous computation of first-passage time PDFs. 

In particular, we shall stress the connections of both 
techniques under the continuum limit approximation. 
This concept — which we will define properly in short — 
can be loosely identified with the limit in which both the 
mean sojourn time fi and the characteristic jump magni- 
tude tend to zero — note that, by contrast, in this paper 
we have considered a situation wherein sojourn times are 
moderate or even large. 
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The major benefit of the continuum assumption is that 
it allows to obtain general results on the basis of limited 
knowledge of the jumping time and size PDFs, even when 
these distributions do not have all their moments well de- 
fined. The major drawback within our set-up is that as 
fi — > the variable r, the time elapsed from the last 
known jump, tends to zero as well, and the distinction 
between T b (x,r) and T b (x) becomes irrelevant. There- 
fore, any comparison between the two methods must be 
focused on how the object T b (x) is obtained. 

Let us assume, for instance, that ip(s) ~ 1 — /is + o(s) 
and h(s) ~ 1 — fc|s| Q +o(|s| a ) when s —J- for certain con- 
stants fi, the mean sojourn time, and k. (Note however 
that it is also possible to consider the continuum limit in 
the case in which the mean sojourn time does not exist. 
We are just giving an illustrative example. For a more 
exhaustive analysis of this topic see H3-) 

To be more explicit, let us consider the case 



h(u) 



k 



(u>0), 



so that, 



h{s) = e- k ^ ~ 1 - kVs + o(i/s). 



(Al) 



(A2) 



It is well known that for the CTRW process S t , Eq. ([I]), 
the propagator p(u, i)du = Pr{u < St ^ u + du} reads in 
the Laplace-Laplace domain 



P(si,s 2 ) 



l-^(aa) 



(A3) 



as si.2 — > 0. The continuum limit is recovered in this 
case by letting /j, — > 0, k — > with k/p — > K finite. 
The process arising after this limit fulfills the following 
fractional diffusion equation 



d t p{u, t) + K DZp(u, t) = 0, 



(A4) 



where fl„ is the Riemann-Liouville fractional operator 



of order i, and whose solution reads 



Kt K 2 t 2 

=e 4 " 



2v 7TU 3 



(A5) 



Let us now define Ub(x,t) as the probability that the 
process X t has never touched b when it started at x < b 
at the initial time, i.e. 

n b (x, t) = Pr{X t , <b,0^t' ^ t\X = x}. (A6) 

In the present case, as St is an increasing positive process 
and v > 0, it can be computed by direct integration of 
p{u,t): 



IL b (x,t) 



b—x—vt 



p(u, t)du = Erfc 



KH 



2+2 



2\]b — x — vt 



Now one can obtain Tb(x) through 



(A7) 



T b (x)= td t [l-IL b {x,t)]dt= IL b (x,t)dt 

Jo Jo 




Alternatively, if we use the direct approach followed in 
this paper our results imply that the mean arrival time 
follows by inversion of the Laplace expression in Eq. (TT2")) 
when /x — > 0, k — > with k/ ' \i — ► K finite 



vs 2 l _ i/j(sv)h(s) vs 2 + Ksi 



(A9) 



that is 



^ rc+ioo e s(b— x) 
T b{x) = — I ; Tds. 



27r * Jc-too VS 2 + KS2 



(A10) 



Upon evaluation of this integral the result (|A8|) is recov- 
ered. 
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